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ABSTRACT 

We calculate the rate of photoevaporation of a circumstellar disk by energetic 
radiation (FUV, 6eV < hu <13.6eV; EUV, 13.6eV < hu <0.1keV; and Xrays, 
hu > O.lkeV) from its central star. We focus on the effects of FUV and X-ray 
photons since EUV photoevaporation has been treated previously, and consider 
central star masses in the range 0.3 — 7M0. Contrary to the EUV photoevap- 
oration scenario, which creates a gap at about ~ 7 (M^/IM©) AU and then 
erodes the outer disk from inside out, we find that FUV photoevaporation pre- 
dominantly removes less bound gas from the outer disk. Heating by FUV photons 
can cause significant erosion of the outer disk where most of the mass is typi- 
cally located. X-rays indirectly increase the mass loss rates (by a factor ~ 2) by 
ionizing the gas, thereby reducing the positive charge on grains and PAHs and 
enhancing FUV-induced grain photoelectric heating. FUV and X-ray photons 
may create a gap in the disk at ~ 10 AU under favourable circumstances. Pho- 
toevaporation timescales for ~ IM© stars are estimated to be ~ 10^ years, 
after the onset of disk irradiation by FUV and X-rays. Disk lifetimes do not vary 
much for stellar masses in the range 0.3 — 3M0. More massive stars (> TMq) 
lose their disks rapidly (in ~ 10^ years) due to their high EUV and FUV fields. 
Disk lifetimes are shorter for shallow surface density distributions and when the 
dust opacity in the disk is reduced by processes such as grain growth or settling. 
The latter suggests that the photoevaporation process may accelerate as the dust 
disk evolves. 
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1. Introduction 

Circumstellar disks are closely associated with star formation, with both strong observa- 
tional evidence for their ubiquitous presence around young stars, and theoretical support as a 
natural consequence of angular momentum conservation in rotating, star- forming cores (e.g., 
Beckwith et al. 1990, Shu et al. 1994). Disk lifetimes are short compared to stellar lifetimes, 
and dust disks around solar mass stars are estimated to typically last < 3 Myrs (e.g., Haisch 
et al. 2001). Gas disk lifetimes are not as well constrained observationally (< 10^ yrs, Zuck- 
erman et al. 1995, Pascucci et al. 2006) but are presumably similar to dust disk lifetimes. 
These timescales set an upper limit on the time available for the assembly of planetary sys- 
tems. Furthermore, dust emission studies that are commonly used to probe disk evolution 
indicate that disks transition from an optically thick stage characterised by strong accretion 
and high ratios of total infrared to bolometric luminosities {Lm/Lboi ~ 10~^, classical T 
Tauri star, CTTS) to a weakly accreting phase with optically thin dust {Ljn/Lboi < 10~^, 
weak-hne T Tauri star, WTTS) (e.g., see recent review on dust disk evolution by Watson et 
al. 2007). Based on the statistics of these two classes of objects, the transition epoch is esti- 
mated to be very brief, ~ O.lMyr, and may involve optically thin inner disks with optically 
thick outer regions (Hillenbrand 2008, Cieza 2008). Optically thick gas and dust disks thus 
originate during the star formation process, undergo significant evolution over a few Myrs to 
perhaps also form planets, and are rapidly destroyed during the planet formation process. 

Various mechanisms have been proposed to date in order to explain the nature of disk 
dispersal (see HoUenbach et al. 2000 and Dullemond et al. 2007 for recent reviews). Disk 
matter builds planetary systems. However, it has long been recognized from considerations 
of the total gas and solid content of the solar system and by comparing initial disk masses 
(~ O.IM^) to the observed inventory of planet mass (~ lO^^M^) around extrasolar systems 
that much of the mass is, in fact, physically removed from the disk either to spiral into 
the central star or to disperse back into the interstellar medium (e.g., Hayashi et al. 1985). 
Viscous accretion is known to occur in the disk and allows transport of mass onto the central 
star, but is insufficient as a dispersal mechanism in itself. Viscosity also spreads the disk, 
as a small amount of matter carries angular momentum outward, and viscous disks do not 
completely disappear at ~ 10 Myr, as observed (gas masses < 10^^ Mj, e.g., Zuckerman 
et al. 1995, Hollenbach et al. 2005, Pascucci et al. 2006; and dust masses < 10^"^ Mj, 
e.g., Beckwith et al. 1990, Haisch et al. 2001). The slow decline in mass of a viscously 
evolving disk (M^ oc t~^, e.g. Hartmann et al. 1998) is moreover inconsistent with the 
observationally inferred rapid transition from the CTTS to WTTS phases. One of the most 
promising mechanisms heretofore proposed for destroying disks has been photoevaporation, 
where the high energy radiation from a young star heats gas at the disk surface and results 
in escaping mass flows from the outer regions of disks. In this scenario, viscous evolution 
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causes the inner disk to accrete onto the central star and to spread, while photoevaporation 
removes the outer disk mass reservoir. 

A number of researchers have worked on the problem of disk photoevaporation in various 
contexts. HoUenbach et al. (1994, hereafter HJLS94), Yorke & Welz (1996), and Richhng 
& Yorke (1997) examined the photoevaporation caused by the extreme ultraviolet (EUV; 
hu > 13.6 eV) photons from the central star. Johnstone et al. (1998), Storzer & Hollenbach 
(1999) and Richling & Yorke (1998, 2000) modeled photoevaporation of disks around low 
mass stars caused by the EUV and the far ultraviolet (FUV; 6 < hu < 13.6eV) fluxes from 
nearby massive stars. Alexander et al. (2004) considered X-rays from the central star and 
found that they were unlikely to lead to signiflcant photoevaporation rates (also see Ercolano 
et al. 2008). Clarke et al. (2001) were the first to combine EUV-induced photoevaporation 
from a central solar mass star with turbulent viscosity to follow disk evolution. They found 
that EUV photoevaporation produces a gap at a characteristic radius ~ 3 — 10 AU once the 
viscous accretion through the disk has declined to a rate M^cc ^ 10~^ M0 yr^^ below the 
photevaporation rate from that region for typical EUV flelds (~ 10^^~^^s~^). The timescale 
for producing the gap was estimated as ~ a few Myrs for sufficiently high EUV flelds and 
small initial disk masses. At this point the inner disk rapidly accretes onto the central star 
on a viscous timescale (~ 0.1 Myr), followed by the rapid photoevaporation of the outer 
disk by the unattenuated EUV flux from the central star (Alexander et al. 2005, 2006 a,b). 
Matsuyama et al. (2003), Font et al. (2004) and Ruden (2004) have also modeled disk 
dispersal by a combination of photoevaporation and viscous spreading and accretion. 

Photoevaporation due to the irradiation of a disk by FUV photons from a central low 
mass star has not been studied earlier. FUV and X-ray luminosities of young low mass stars 
are signiflcant and better determined than their EUV luminosities, which cannot be directly 
observed. In young low mass stars, FUV, EUV and X-ray fluxes are far in excess of older 
main sequence low mass stars like the Sun. This is due to both enhanced magnetic activity 
on stellar surfaces, which leads to much higher chromospheric activity and production of 
energetic photons, and to the presence of substantial accretion (in accretion columns, Gull- 
bring et al. 1998, Calvet & Gullbring 1998) onto the central star from the disk, which leads 
to accretion shocks producing energetic photons. 

Although FUV, EUV, and X-rays are produced near the stellar surface, they may be 
absorbed before reaching the disk surface. Alexander et al. (2005) have shown that the EUV 
produced in the accretion shock is not likely to penetrate the accretion column. Disk winds 
accompany accretion with a wind mass loss rate that scales as ~ xO.l the accretion rate. In 
order for the chromospheric EUV to penetrate the disk wind, the accretion rate needs to be 
less than ~ 10~^ Mq yr~^ (see Hollenbach & Gorti 2008). FUV and X-ray photons require 
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much larger column densities than the EUV to be absorbed and therefore begin to penetrate 
the protostellar wind and photoevaporate the disk once accretion rates are < 10~^ Mq yr~^ 
(HoUenbach & Gorti 2008), long before the EUV photons from the stellar chromosphere can 
impact the disk surface. 

In this paper, we include the effects of FUV and X-ray radiation from the central star on 
disk photoevaporation. This is a preliminary study and our main goal is to assess the impor- 
tance of FUV radiation on disk evolution and lifetimes. Wc therefore make many necessary 
simplifications, notably the use of steady-state, static models of the chemical abundances, 
gas and dust temperatures, and density structure followed by a simple "streamline" analysis 
of the photoevaporative flow. We investigate the relative importance of optical, FUV, EUV 
and X-ray radiation from the central star in inducing photoevaporation by using our recently 
developed disk models (Gorti & HoUenbach 2008, hereafter GH08). We present models of 
disks around central star masses of 0.3 — 7 M© that span the range of current observational 
disk detections (e.g., Luhman et al. 2005, Fuente et al. 2006, Manoj et al. 2007). We 
consider disks in their initial stages, with moderate grain growth and power law surface den- 
sity profiles. We calculate photoevaporative mass loss rates from the disk as a function of 
relevant model parameters such as the dust opacity per H nucleus, the FUV luminosity, the 
X-ray luminosity and the power law exponent of the surface density distribution. Although 
our results are for a snapshot in time, we quahtatively estimate the effects of viscosity and 
discuss disk evolution and dispersal. The paper is organized as follows. In §2 we present a 
brief overview of disk evolution prior to the onset of photoevaporation and describe our disk 
models and photoevaporation calculations. In §3 we discuss our results and their implica- 
tions for disk lifetimes, and describe disk evolution as a function of various input parameters. 
We summarize and conclude in 54. 



2. Model Description 

2.1. Onset of Photoevaporation 

We consider the stage in the star formation process when disk illumination by FUV 
photons and X-rays is likely to commence. During the gravitational collapse of the proto- 
stellar core and the formation of a central star and circumstellar disk, accretion of mass onto 
the star proceeds mainly through the disk. In this earliest phase of stellar evolution, average 
accretion rates are very high (> 10~^ Mq/yv) and there is complete obscuration of the UV 
and X-rays from the central star by infall onto the inner disk and the star. During this 
phase, disk accretion and transport of angular momentum may be accomplished by global 
gravitational phenomena such as spiral density waves, as the disk hovers on the brink of 
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gravitational instability. It is likely that accretion rates through the disk and onto the star 
are comparable with infall rates onto the disk. As infall onto the disk ceases, and the disk 
continues to accrete onto the star, the ratio of disk mass/stellar mass decreases and the 
self-gravity of the disk becomes less important. Turbulent viscosity may then provide the 
means of angular momentum transport and the disk mass declines, the disk spreads and the 
accretion rate on to the star steadily drops (Najita et al. 2007, Bouvier et al. 2007, Watson 
et al. 2007). Disk illumination by FUV and X-rays commences once they can penetrate the 
protostellar wind originating in the inner disk. 

Accretion onto the star and stellar magnetic activity provide sources of energetic pho- 
tons. There is a significant output of energetic photons from the star due to the accretion 
hotspots on its surface (FUV, and perhaps EUV and soft X-rays, e.g. Calvet & GuUbring 
1998, Robrade & Schmitt 2006, Giidel et al. 2007) and due to chromospheric activity (bulk 
of X-ray emission, e.g, Preibisch et al. 2005 and likely EUV, Alexander et al. 2005). 

Accretion is also accompanied by protostellar winds originating from the inner disk that 
are proportional to the rate of accretion (M^ ~ 0.05 — O.lMacc, White & Hillenbrand 2004). 
Considerable amounts of high energy radiation from the star are initially absorbed at the 
base of the accretion column and in the wind (e.g, MuzeroUe et al. 2001). As accretion 
diminishes, the opacity of the stellar outflow and the wind from the inner disk is low enough 
that the outer regions of the disk arc irradiated by stellar optical, FUV and X-ray photons. 
Accretion shock-generated FUV and soft X-ray photons are observationally inferred to at 
least partially penetrate the accretion column (e.g., Muzerolle et al. 2003, Robrade & Schmitt 
2006), and have been measured around many young stars (e.g., lUE observations for FUV, 
Valenti et al. 2003). The hydrogen nucleus column density of gas necessary for absorption 
of X-rays (Nh ~ lO^^cm"^) and FUV (Nh > lO^^cm"^, dependent on dust opacity in the 
protostellar wind) is larger than for EUV absorption (Nh ~ lO^'^cm"^). X-rays and FUV 
photons therefore penetrate the wind and begin to irradiate the disk at earlier epochs than 
EUV photons. Hollenbach & Gorti (2008) show that the FUV and X-ray photons begin to 
penetrate the wind once the accretion rate drops below ~ 10^^ M0 yr^^(>; 0.1 — 1 Myr), while 
EUV photons penetrate the wind only when the accretion rate drops below ~ 10~^ M0/yr 
(> l-3Myr). 



2.2. Model Assumptions and Input Parameters 

The main input parameters to the model arc the stellar mass, radiation fleld, the gas 
phase abundance of elements, the surface density distribution of gas and dust, and the dust 
properties in the disk. In this paper, we calculate the disk structure and photoevaporative 
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mass loss rates at a characteristic time (~ 1 Myr) during the planet-forming disk evolutionary 
epoch. From the mass loss rates, we estimate disk dispersal times. We therefore use input 
parameters such as stellar radiation field and dust properties that are representative of this 
characteristic time. We list the assumed properties for a fiducial IMq star case in Table [T] 
and for different stellar masses in Table [2l 

Stellar Parameters We assume stellar parameters such as radius, effective temperature 
and bolometric luminosity as appropriate for a pre-main-sequence star of age ~ 1 Myr and 
use the evolutionary tracks of Siess et al. (2000). We consider a range of stellar masses 
from 0.3 — 7 Mq. We do not consider lower mass central objects, although infrared excesses 
have been noted around objects with masses as low as O.IMq. It is rather uncertain as 
to whether stars more massive than ~ TMq harbor disks around them (e.g., Fuente et al. 
2007, Manoj et al. 2007). In fact, we will show that massive stars (> 7M0) photoevaporate 
their disks on very short timescales (~ 10^ years). The X-ray luminosities of stars are fairly 
well determined and we use observational data to guide our choice of X-ray luminosity as 
a function of mass {Lx ~ 2.3 x 10^°{M^/MQy-^^eTg s"^ for < 3Mq and Lx ~ lO'^L, 
for > 3Mq; Flaccomio et al. 2003, Preibisch et al. 2005). The EUV luminosities of 
stars are very poorly known and at best, indirectly determined (e.g., Bouret & Catala 1998, 
Alexander et al. 2005). Here the EUV luminosity (that is not from the photosphere) is 
assumed to be chromospheric in origin like the X-rays, and of a similar strength and scaling 
with mass. FUV luminosities of stars are well-studied observationally (e.g., lUE, Valenti 
et al. 2003, FUSE, Bergin et al. 2003, Herczeg et al. 2004) and theoretically believed 
to arise due to both activity in the chromosphere and accretion hotspots on the surface of 
the star. We consider both components for our FUV spectrum. We calculate accretion- 
generated FUV luminosities from mass-correlated accretion rates (e.g., Gullbring et al. 1998 
for accretion generated FUV, MuzeroUe et al. 2003 for mass-dependent accretion rates; see 
GH08 for details of our procedure). The chromospheric component to the FUV flux has 
an FUV luminosity given by logLpjjy/L^ = —3.3 (data from Valenti et al. 2003), similar 
to the scaling for X-ray luminosity. However, as young stars are known to be variable and 
as disk accretion rates decline with time, we treat the FUV and X-ray luminosities as free 
parameters and vary them over several orders of magnitude, keeping other parameters fixed. 

Disk and Dust properties We assume the initial disk mass to be proportional to the 
mass of the star with M^isk ~ 0.03M^, (e.g., Andrews & Williams 2005). Our disk extends 
from 0.5 to 200 AU and the surface density distribution is generally assumed to follow a 
power law, S(r) oc r~^. In one case we assume a power law surface density distribution 
to probe the sensitivity of the results to the power law exponent. 
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Dust in the disk is assumed to have marginally evolved at epochs of consideration, 
and we assume a power-law size distribution for the dust with a range from 50A— 20/im, 
corresponding to a factor of 10 reduction in opacity compared to interstellar dust, if the 
gas/dust mass ratio in the disk is assumed to be 100. We also consider further reductions in 
opacity, which we treat by retaining the same size distribution but increasing the gas/dust 
mass ratio. Our modeled dust has a mixture of chemical compositions (GH08). Therefore, 
at the disk surface where the dust is heated by stellar photons, the temperature of a dust 
grain is a function of both its size and composition. 

Disk Thermo-chemical Models Photoevaporative mass loss rates depend sensitively on 
the density and temperature of the gas in the flow, and in order to calculate X-ray and FUV- 
initiated photoevaporation rates, detailed thermo-chemical disk models are required. EUV 
photoevaporative flows, on the other hand, are isothermal due to thermal balance processes 
in ionized gas which set the gas temperature to be nearly constant at ~ lO^K. Heating and 
cooling of the predominantly neutral gas below the ionized surface layer is complex, with dust 
collisions. X-rays, FUV, photo-reactions and chemistry all being important in determining 
the density and temperature structure of the gas. We use our recently developed thermo- 
chemical, steady-state, disk models to determine the gas density and temperature structure 
in the disk (GH08). 

We briefly summarize our disk model here and refer to GH08 for more details. Our disk 
models solve for chemistry, thermal balance and impose vertical hydrostatic equilibrium (all 
self-consistently) to separately calculate the density and gas and dust temperatures as a 
function of spatial location in the disk. We consider heating of the gas due to X-rays, grain 
photoelectric heating by PAHs and small grains, cosmic rays, exothermic chemical reac- 
tions, formation heating of H2, collisional de-excitation of vibrationally excited (by FUV) 
H2, photoionization of carbon, and collisions with warmer dust grains. Cooling of gas is 
by line emission from atoms, ions and molecules, and by collisions with cooler dust grains. 
Our chemical network is moderate and focused towards including species that are dominant 
coolants in the disk (enabling an accurate determination of the gas temperature). We con- 
sider 84 species (ions, atoms and molecules) of the elements H, He, C, O, Ne, S, Mg, Fe, Si 
and Ar and ~ 600 chemical reactions, including photoreactions and ionizations caused by 
cosmic rays. Our models neglect ice formation on grain^ and calculate steady-state chem- 
istry. We treat gas and dust independently, and allow for different spatial distributions, 
though we generally consider them well-mixed throughout the disk. Dust radiative transfer 



^In GH08 we show that in the surface layers where photoevaporative flows originate, a combination of 
thermal and photodesorption prevent substantial accumulation of ice on grains. 
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for optically thick dust disks is simplified, so as to keep the numerical disk model compu- 
tations tractable. We use the two-layer model for dust with some modifications (Chiang & 
Goldrcich 1997, Dullemond, Dominik & Natta 2001, Rafikov and de Colle 2006). We include 
the effects of background infrared radiation due to dust on the level populations of the gas 
species (see Hollenbach et al. 1991). We then calculate the gas temperature, density and 
chemical structure as a function of spatial location (r, z) throughout the disk. 

Assumption of steady-state evolution We emphasize that, in this paper, our entire 
treatment is steady-state. We use stellar properties and radiation fields and disk surface 
densities at a characteristic age (arbitrarily set at ~ 1 Myr), although these are expected to 
evolve with time as the star evolves. Our choice for the age of the system is motivated by the 
observational evidence for disk hfetimes of a few Myrs (e.g., Haisch et al. 2001, Silverstone 
et al. 2006, Siciha-Aguilar et al. 2006) and the expectation that it might be a few 10^ years 
before the star is visible to the disk surface because of the heavy accretion of matter from the 
dense molecular core onto the star and inner disk (Shu et al. 2000). The characteristic age of 
~ IMyr therefore approximately corresponds to the onset of FUV and X-ray irradiation of 
the disk, following penetration of the massive protostellar wind arising from the inner disk. 
Disk properties also evolve, with increased dust settling and grain growth as time progresses, 
but in these static models we fix our dust properties, such as the opacity in the disk and 
the gas/dust mass ratio. The surface density distribution also evolves under the combined 
effects of photoevaporation and viscous evolution (e.g. Alexander et al. 2006). However, 
viscous evolution, if dominant and characterized by a constant viscosity parameter a, will 
tend to maintain a surface density distribution that is approximately proportional to r~^. 
Our aim in this initial study is to estimate the importance of FUV-induced photoevaporation 
in disks and in this preliminary work we ignore all time- dependent processes. Instead, we 
focus on deriving average rates of photoevaporation, by using stellar and disk properties at 
an intermediate epoch in disk evolution. We also use an average EUV field through the 
disk evolution, although EUV photons may not illuminate the disk until later epochs when 
the mass loss rate has dropped to < 10~^ Mq/jt. In a future paper, we will couple viscous 
evolution and photoevaporation in a time-dependent calculation and use evolving stellar 
fluxes (Gorti, Dullemond & Hollenbach 2008, in preparation). 

We refer the reader to GH08 for a more thorough discussion on all the above assumptions 
and choice of input parameters. 
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2.3. Photoevaporation Theory 



Photoevaporativc flows are driven by thermal pressure gradients in the gas, and therefore 
the gas temperature and density distribution (obtained from our disk models) determine the 
mass loss rates. Although our models are steady-state and the disk evolves with time, we 
assume the steady state temperatures and densities do not deviate significantly from the 
time-dependent values of these physical variables at the surface where the flow originates. 
Photoevaporativc flows from FUV and X-ray heated surfaces are launched subsonically from 
regions high in the disk which suffer little extinction to the central star, as will be shown 
in §3. Therefore, heating and chemical timescales are short, of order the timescales for UV 
photodissociation and X-ray ionization (~ 10^ years at 100 AU), compared to dynamical 
timescales in the flow {tdyn ~ r/vjiow ~ ''"/(O.lCs) ~ 10^^^ years at 100 AU). Therefore, the 
assumptions of steady-state are approximately valid and our calculations provide reasonable 
estimates of the mass loss rates. The hydrodynamics assumes an isothermal flow from the 
launch point to the sonic radius, and spherical geometry. Note that the photoevaporativc 
flows we envisage from the surfaces of disks are distinct from the usually bipolar protostellar 
winds present during the early stages of star formation. We do not treat any possible 
interaction of the protostellar and photoevaporativc winds, which is a complex phenomenon 
beyond the scope of the present work. However, we expect the interaction to result in 
shearing layers which may possibly enhance the photoevaporativc mass loss rates from the 
disk (Matsuyama, Johnstone & HoUenbach, in preparation). 

The rate of photoevaporation (at a given radius r) is dominated by the layer at a height 
Zfiow where the combination of density and temperature yield the maximum mass efflux. 
From our calculation of the density and temperature structure of the gas due to heating by 
FUV, EUV, X-rays and other processes, we flnd the height Zfiow where this maximum occurs 
and thereby obtain the photoevaporation rate as a function of disk radius. As shown by 
HJLS94, a characteristic gravitational radius enters all discussions of photoevaporation; 
rg is the radius where the sound speed or the thermal speed of hydrogen is equal to the 
escape speed from the rotating disk. For ionized gas, Vgji ~ 7(1O^K/T)(M*/1M0) AU. In 
predominantly neutral gas heated by FUV and X-rays 



Note that the gas temperature T varies with z, so that it is possible at a flxed r to have 
certain heights z where r < Vg and some z where r > Vg. If r > r^, then the potential mass 
flux from a given height z is given by 




(1) 



E(r) 



rfS(r) 



// n{r,z) Cs{r,z) 



(2) 



dt 
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where n is the mean particle mass per H nucleus, n is the hydrogen nucleus density, and Cg 
is the sound speed (see HJLS94). 

For lower gas temperatures so that r < Vg, Adams et al. (2004) showed that subsonic 
flows can be launched. In fact, substantial mass fluxes occur for r > where ~ 0.1— 0.2rg 
(also see Begelman et al. 1993, Liffman 2003, Font et al. 2004). The flow accelerates through 
the sonic point (where the sound speed is attained) and the photoevaporative flow is launched 
with the mass flux given by (see Adams et al. 2004) 

fmsCs{r, z) [ — ) (3) 



at \ r 

where Us the density at the sonic radius and is the sonic radius. If we make the assumption 
that the gas is isothermal from the launch point to the sonic radius, and that r,, ^ r so that 
spherical symmetry is a reasonable approximation, then 



= n(r,z)exp(-^(l--)^) (4) 



and Tg is given by 



1/2^ 



9 . 

for r < r^/S. A Parker wind solution has = rg/2. For r > rg/2, the flow rapidly goes 
through a sonic point near the base and = r. For r^/S < r < rg/2 we linearly extrapolate 
between = Vg/A at r = r^/S and = rg/2 at r = rg/2. Using Eqs. ([2HS]) and the density 
and temperature distribution from our disk models, the mass loss rate at each radius is 
determined by the z layer with the highest E, which could be heated by FUV, EUV, X-rays 
or any of the other heating processes (§2.4). 

Note that the above analysis assumes that the flow is isothermal from r to r^, whereas 
a more accurate hydrodynamical analysis may produce a variation in temperature along a 
flow streamline. We postpone a more complete analysis to future work and at present check 
the validity of our solutions by following a radial ray from the star through the disk from 
the flow surface to the sonic radius. In fact, we find from our thermal balance calculations 
that the temperature is likely to increase marginally along the flow streamline, increasing 
the photoevaporation rate from that given by Eqs. 



In the case of EUV-driven photoevaporation, because of the isothermal nature of the 
ionized region, the flow originates just above the ionization front where the density is greatest. 
At this height -Z/zow, the density is given by (see HJLS94) 
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where = 2.53 x 10 ^^cm^s ^ is the case B recombination coefficient of hydrogen, rgjj = 
jj is the gravitational radius, (pEuv s ^ is the EUV photon luminosity of the star 
and the exponent p is equal to 1.5 for r < Vg^j and 2.5 for r > Vg jj. The mass loss rate (or 
equivalently r^S) peaks at ~ // and decreases beyond // as the density at the base riji 
decreases rapidly with r (HJLS94). 

The photoevaporation of gas occurs at the surface of the disk and carries with it small 
dust particles. In the absence of turbulence, where dust is allowed to freely settle, the 
criteria for dust removal is that the terminal speed of a settling dust particle must be less 
than the speed of the downward moving photoevaporative front. Such a calculation shows 
that even micron-sized dust particles will settle fast enough to avoid flowing outward with 
the evaporating surface gas. In other words, the dust abundance in the upper surface layers 
is extremely low. In the other extreme, one can assume that turbulence mixes the dust 
vertically so that no settling occurs. In this extreme, the criterion for dust removal is that 
the ram pressure of the gas moving upward past the dust is sufficient to accelerate the dust 
to escape speed. This criterion results in dust with sizes < mm being coupled to the gas 
and evaporating while larger dust particles or objects remain behind in the disk (Adams 
et al. 2004). Since larger particles are decoupled to the flow and may remain in the disk, 
photoevaporation can decrease the ratio of gas surface density to the dust surface density 
in disks (see Throop & Bally 2005 and Youdin & Shu 2002 for a discussion of the possible 
consequences) . 



3. Results: Photoevaporation and Disk Lifetimes 

3.1. Standard Disk Model for a IMq Central Star 

We ffist describe in detail the results for our standard fiducial disk around a IM© star. 
(We refer the reader to GH08 for a discussion of the disk structure and chemistry.) Figured] 
shows the position of the hot dust surface layer (Ay = 1 to the star at z = zcg in the 
two-layer dust model, e.g., Chiang & Goldreich 1997), the ionization front (where the EUV- 
induced flow would originate) and the FUV/X-ray-induced photoevaporative flow surface for 
the entire disk. Although we discuss the ionized flow surface and the FUV and X-ray heated 
neutral flow surface as being distinct for clarity, the true location of the base of the mass flow 
is at the height where the maximum mass loss (due to EUV, or due to FUV and X-rays) 
occurs. At the very surface (Ay < 10^^) of the entire disk, EUV heating (see HJLS94) 
proceeds via the photoionization of hydrogen and leads to a nearly isothermal 10^ K H II 
region from the ionization front to the surface. The gas rapidly becomes neutral and colder 
below the ionization front. The EUV photoevaporation rate therefore peaks at the base of 
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the EUV-heated layer, where the density is highest. The neutral flow surface (z/zow), where 
heating is by FUV/X-rays, is seen to generally trace zcgi but lies at slightly higher z where 
Av ~ 0.01 for r < 10 AU, Ay ~ 0.1 - 0.2 for r ~ 10 - 60 AU, and Ay ~ 0.3 beyond 60 AU. 
This behaviour can be qualitatively understood as follows. The mass loss rate increases with 
increasing temperature (T) and density (n), and whereas T increases with z due to greater 
penetration of FUV and X-ray photons, n decreases with z to maintain vertical hydrostatic 
equilibrium. The mass loss rate therefore naturally peaks at a height in the gas temperature 
distribution where both T and n are relatively high. (See GH08 for a discussion of the 
vertical density and temperature structure in the disk). This approximately corresponds to 
the height zcg where the optical depth due to dust is of order unity and where FUV photons 
from the star begin to be attenuated, i.e. when Ay ~ 0.1 — 0.5. 

The gas temperature and density at Zfiowi which are the main parameters that disk mass 
loss rates depend upon, can be approximately estimated by simple considerations. In order 
for matter to escape the disk, the gas temperature at the base of the flow has to necessarily 
satisfy the criterion r > r^. with Tcr ~ 0.1 — 0.2G'M*/Cg (see Eq. |5l also Liffman 2003, 
Adams et al. 2004). This implies a criterion for a minimum temperature at zjiow, T{r) > 
15OO(lOy4f//r)(M*/lM0) K. As the mass loss rate is also a function of gas density which 
decreases with z, the maximum mass loss typically occurs at the lowest height where the 
minimum flow temperature criterion can be satisfied and the density is still high. Thus, the 
temperature at Zfiow is approximately given by the relation T(r) > 15OO(lOAf//r)(M=K/lM0) 
K. Simple considerations show that because the flow surface is situated close to Ay ~ 0.3 for 
most of the disk, the density at Zfiow lies in the range ~ 10^ — lO^cm"'^. The hydrogen column 
density from the star to the Ay ~ 0.3 surface is ~ 6 x lO^^cm"^ for our standard model. 
The local density is ~ 6 x 10^^cm~^/r ~ 4 x 10^(100AU/r)cm~^. Figure [1] demonstrates 
the results of a more detailed analysis using our disk models and shows the density and 
temperature of the gas at the FUV/X-ray flow surface. The gas temperature where the 
mass loss originates declines from ~ 3000K in the inner (< lAU) disk to ~ lOOOK at 10 
AU and ~ 80K at 100 AU in the outer disk. The density at the base of the FUV/X-ray 
heated flow does not change appreciably, and remains ~ 10® — 10'' cm"'^ for disk radii 10 
AU< r <200 AU. 

Heating at the neutral flow surface in the inner disk (~ 3 — 30 AU) is mainly due to 
FUV photons (~ 70%), and there is some contribution from chemical processes such as H2 
formation heating and X-rays. Cooling is due to [OI]63/im line emission, H2 vibrational 
and rotational lines, and dust collisions. In the outer disk (r > 30AU), FUV heating still 
dominates (~ 80%) and there is some X-ray heating. Cooling at the flow surface is mainly 
due to [01]. Gas temperatures at the flow surface in the outer disk can be below the mean 
dust temperature, but in general, the gas is nearly always warmer than the dust at the flow 
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surface so that collisions of gas and dust cool the gas. 

Having determined the location of the flow, Zfiow, at every radius as the height with 
maximum S(r), (Eq. [3]), we next estimate the mass loss rate from the disk. Figure [2] shows 
the mass loss rate Mpe for our standard disk model (Model S, solid line) as a function of disk 
radius, for the resulting net (neutral or ionized) flow. Mpe = 27rr^S(r) is approximately the 
mass loss rate between r/2 and 3r/2, i.e, Ar ~ r centered on r. In the inner disk (r <3 AU), 
mass loss rates driven by EUV radiation are higher whereas in the outer regions of the disk 
(r > 3 AU) FUV/X-ray photoevaporation is found to dominate. At r ~ 3 AU, there is a rapid 
rise in Mpe as Zfiow shifts from the ionized region to the neutral, dense FUV and X-ray heated 
region. Although mass loss for EUV photoevaporation peaks at ~ rgjj ~ 7 AU, FUV/X- 
ray induced mass loss rapidly begins to dominate once r > 3 AU for our standard FUV 
luminosity which corresponds to an accretion rate of 3 x lO~^M0yr~^ . Gas here is heated 
by FUV and X-rays to temperatures of ~ 2000K, and the higher densities in this region 
result in higher photoevaporative rates than in the lO^K, lower density gas in the ionized 
region. There is another shift in the location of Zfiow at ~ 10 AU, where the flow surface 
shifts from Ay ~ 0.01, FUV (80%) and X-ray (20%) heated gas to denser, predominantly 
FUV-heated gas at Ay ~ 0.1. At this radius. Fig. [2] shows a corresponding kink in Mpe. 
Beyond ~ 10 AU, Zfiow remains approximately at an Ay of 0.3, as discussed above. However, 
the surface gas temperature rapidly falls with radius as the FUV/X-ray flux decreases away 
from the star and Mpe drops sharply at ~ 25 — 40 AU. The exponential term in E, e~^''^'^^ 
(see Eq. 12]), tends to dominate the mass loss rate in this region. Since Vg/r oc (Tr)~^, and 
because T falls faster than (see Fig. [T]), the mass loss rate at ~ 25 — 40AU decreases 
with increasing r. Beyond ~ 40 AU, the gas temperature at the flow surface is a decreasing 
shallow function of radius varying as ~ r-o-^ and the binding energy of the disk gas begins 
to decrease, permitting cooler gas in these regions to escape in neutral flows (decreasing 
Tg/r factor in Eq. [3]). Here, although the e"^"/^^ term and E are relatively constant with 
r, Mpe oc r^E and therefore increases with r. There is a small kink in the location of Zfiow 
at ~ 60 AU(Fig. [1]), where the peak in E shifts from Ay ~ 0.18 to Ay ~ 0.34 but E and 
hence Mpe varies only slightly with z at each radial zone in this region, and Mpe is relatively 
smooth with r. 

The peak mass loss rates occur at the outer edge of the disk and are induced by FUV 
photoevaporation, which therefore sets the disk lifetimes. Typical surface density distribu- 
tions in disks are shallow (E(r) oc r~", a ~ 0.5 — 1) such that most of the disk mass is at larger 
radii (e.g., Andrews & Williams 2007). FUV photoevaporation is therefore responsible for 
the depletion of the disk mass reservoir. A simple estimate can be made of the characteristic 
timescale tpe for FUV-induced photevaporation, by considering the flow parameters in the 
outer disk. As shown earlier, zjiow ~ ^cg- For our assumed dust opacities, the density at 
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the base of the flow, n{r, Zfiow), can be approximated as ?7.(100AU, Zfiow) = 4 x 10^ cm~^. At 
these densities, thermal balance and chemistry calculations determine the gas temperature 
as ~ 80K. Equation [3] then obtains S ~ 2 x 10~^^ gm cm~^ s~^. For a 200 AU disk, with mass 
0.03 Mq and S ~ the surface density at 100 AU is ~ 1 gm cm~^ s~^. The characteristic 
timescale for FUV-induced photoevaporation at 100 AU, where most of the mass resides, is 
therefore short, with tpe = S/E ~ 10^ years. 

While the disk on the whole begins to lose mass as it is irradiated by FUV and X-ray 
photons, the morphological evolution of the disk depends on the relative and local rates 
of photoevaporation and viscous evolution. Viscosity also plays an important role as it 
simultaneously spreads the photoevaporating disk. However, viscous timescales in the outer 
disk are typically longer than the FUV photoevaporation timescales, and viscous evolution 
cannot replenish the photoevaporating disk at large radii. Therefore, the mass reservoir 
here is rapidly drained. In Figure [3l we show the characteristic photoevaporation timescale, 
tpe and the viscous timescale for comparison. tpe{r) is calculated as S(r)/S(r), and 
viscosity replenishes mass locally on the viscous timescale ~ 2r^/3z/, where the viscosity 
u = aCsH{r), a is the viscosity parameter, Cg is the local sound speed and H{r) is the disk 
scale height. We use our disk structure calculations to obtain H{r) and Cs- Our initial 
assumptions on the surface density distribution with radius and steady mass accretion rate 
{Mace = 3vrz/S) in the disk together specify the value of the viscous parameter, a (e.g., 
Pringle 1981, Hartmann et al. 1998), which for our fiducial disk is calculated as ~ 0.02. 
Disks may have a range of a, depending on the extent of MRl-induced ionization. In §3.2 
we consider models with the same initial disk mass but varying mass accretion rates (and 
hence FUV luminosities) which correspond to different values of a. Note that the value of 
a does not directly affect Mpe (which depends on the incident FUV luminosity) but does 
determine the viscous timescale t^. We thus calculate t^, and the two timescales are shown 
in Fig. [31 tpe decreases with r in the outer disk, while ti, increases with r. A truncation 
radius can be defined as the radius where these two timescales are equal in the outer 
disk. For r > rt, Mpe is high, tpe < t^, and viscous timescales are too long to replace 
mass loss by photoevaporation. The disk is therefore rapidly truncated to this radius, and 
continues to photoevaporate and viscously evolve while roughly maintaining this size. We 
estimate disk lifetimes by comparing the timescales tpe and t^ (Figure [3]). tpe = t^, at the 
truncation radius rt. Mass loss is faster than viscous spreading for r > in the disk, and 
viscous accretion removes disk mass at r < r^. The disk loses mass over this characteristic 
timescale, tpe = t^ = Tdisk, which we define as the disk lifetime. For our fidicual disk, with 
a constant a = 0.02, we obtain a disk lifetime of 10^ years after the onset of FUV/X-ray 
photoevaporation. 

Gaps can be created at a given radius when tpe{r) < t,y{r), or equivalently, when 
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Mpe{r) > Mace- When tpe > t^, viscous evolution rapidly transports mass from the outer 
disk into the local radial annulus, and as long as the outer disk remains massive, photo- 
evaporation does not create gaps in the disk. For EUV photoevaporation which peaks at 
r ~ rgji (Fig. [2]), when accretion rates in the disk decline to low enough values, a gap 
forms at r ~ Vgji and thereafter the outer disk evaporates from inside-out (Clarke et al. 
2001, Alexander et al. 2006). FUV/X-ray heated flows may also create gaps at ~ 3 — 30 
AU where there is a local peak in Mp^ (also see Fig. [3]). However, for our fiducial disk. 
Mace = 3 X 10^® MqYt~'^ and is somewhat greater than the peak value of Mpg at ~ 12AU, 
indicating sufficient viscous replenishment of mass to prevent gap formation. 

Creating a gap in the inner disk and maintaining the gap against the smearing effects 
of viscosity may be a time- dependent phenomenon for FUV-induced photoevaporation. The 
FUV flux in our standard model (S) is largely accretion generated, but as the disk evolves, 
later lower accretion rate epochs would result in weaker incident FUV flelds, and hence also 
lower Mpe- Maintaining a gap once it is created therefore depends on how tpe changes as 
accretion activity decreases with time. If at any instance of time, tpe < ti, and a gap forms 
at r ~ lOAU, then the gas interior to the gap drains onto the star on the viscous timescale 
here, ~ 10^ years. However the depletion of the inner disk will result in a lower Mace onto 
the star, hence lowering the accretion-generated FUV flux and the resulting Mp^. In the 
absence of accretion, a lower limit to Mpe is set by the chromospherically generated FUV 
flux. If Maec in the region beyond the gap is higher than this value, then the draining of the 
inner disk after the flrst gap creation event would lead to a rapid accretion of the outer disk 
to flll the inner hole created. The star may perhaps undergo periodic "bursts" of viscous 
accretion and rapid photoevaporation from the ~ 10 AU region at these epochs, until the 
accretion rate drops below the lower limit to Mpe that corresponds to chromospheric FUV 
and X-ray photoevaporation, or until perhaps EUV photoevaporation becomes signiflcant. 
Figure [2] shows the mass loss rate for a model (Model FC; labels F and C for FUV and 
chromosphere respectively) where the accretion rate is zero, i.e., the FUV from the star is 
generated entirely by chromospheric activity. In this case, EUV photoevaporation dominates 
the inner disk evolution, up to a large disk radius of almost 100 AU, and presumably would 
maintain (or create) a gap in the disk at ~ 7 AU when Mace ^ 5 x 10~^°MQyr~^. In the outer 
disk beyond ~ 100 AU, heating at the flow surface is by both FUV (50%) and dust collisions 
(50%), and the mass loss rate in the outer disk is still reasonably high, ~ lO~®M0yr^^. 
We show in §3.3 that increasing the X-ray luminosity to >10 times its standard value may 
increase Mpe from FUV/X-rays to a value where it may exceed Mace and hence create a gap. 

Based on the above results, we propose that disk dispersal around a IM© star proceeds 
as follows. In the initial stages of evolution, as the accretion rate and hence the column 
density through the protostellar wind from star to disk surface decrease with time, the 
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FUV photons (with absorption columns ~ 10^^ cm~^) begin to irradiate the disk before 
the wind ceases to be opaque to EUV photons (absorption column ~ 10^° cm~^). In the 
outer regions of the disk, where M^e > Macci photoevaporation dominates and removes disk 
material faster than it can be replenished by viscous spreading. The outer disk regions 
are eroded rapidly and the disk shrinks. For plausible surface density distributions, most 
of the mass of the disk is in the outer regions and much of the disk mass is therefore 
dispersed via FUV-induced photoevaporation. The disk is expected to rapidly shrink to a 
truncation radius (~ 150 AU). Thereafter, viscosity continuously spreads the disk as it 
loses mass due to FUV photoevaporation, approximately maintaining this outer disk size. 
Under favourable conditions (for example, high X-ray luminosities), FUV and X-ray photons 
may create a gap at 3 — 30 AU, and cause periodic episodes of inner disk ablation and viscous 
depletion. Initially EUV photons have very little effect on the disk evolution, but at later 
stages penetrate the wind and begin heating the surface. Finally as the accretion rate 
declines, a gap forms in the disk where Mpe due to EUV is a maximum, at r ~ Tgji. At this 
stage of evolution, most of the disk mass has already been lost from the disk by FUV/X- 
ray photoevaporation. Once the gap forms, the inner disk rapidly disappears on a viscous 
timescale of < 10^ years. EUV, FUV and X-ray photons now directly impact the newly 
formed rim and the outer disk is subsequently rapidly eroded (> 10^ yrs) from inside out 
(as discussed in Alexander et al. 2006). 



3.2. FUV Luminosity 

Gas temperatures at the flow surface of the disk are mainly set by the heating due 
to the grain photoelectric effect initiated by the incident stellar FUV flux. Uncertainties 
and intrinsic variabihty of the stellar FUV field can therefore affect the photoevaporation 
rate from a disk. We investigate the sensitivity of Mpe and Tdisk to our assumed FUV field 
by considering three additional cases. In the earliest stages of evolution, the star actively 
accretes and can have accretion rates much higher than that assumed in our standard model. 
We increase our adopted FUV luminosity by a factor of 10 (corresponding to an accretion 
rate Mace ~ 3 X lO^^M© yr^^) and calculate the disk structure in this case (Model FIO). 
We also consider a case with the FUV luminosity lowered by a factor of 10 (Model FO.l). 
We next consider a disk where accretion has ceased (Model FC). FUV radiation now stems 
only from stellar chromospheric activity and scales with the X-ray luminosity (Wood et al. 
2005). Using archival lUE data (Valenti et al. 2003) we find that the median fractional UV 
excess luminosity for stars with low accretion rates (we use a procedure similar to Calvet et 
al. 2004) is given by ~ Log \juv,excess/'^* = —3.3 and we use this value. For comparison, we 
also consider a case where there is no UV or X-rays and the disk is heated only by stellar 
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optical photons (Model O). All of these models have the same initial disk mass, and hence 
the derived viscous parameter a for each disk scales directly with the mass accretion rate, 
with the fiducial value being 0.02 (Model S). 

Photoevaporation rates (and hence derived disk lifetimes) are quite sensitive to the 
assumed FUV flux and higher (lower lifetimes) for higher FUV fields. Figure H] shows the 
mass loss rate and the corresponding accretion rate for our different runs. At earlier stages 
of disk evolution, high accretion rates and hence high FUV fluxes incident on the disk will 
cause rapid mass loss from the outer regions, depleting disk mass significantly. As the disk 
evolves, both Mace and Mpe decrease, and for the range of FUV luminosities considered here, 
we find that Mpe decreases faster than Mace for accretion-generated FUV. From Model FIO 
to Model S, Mpe decreases by orders of magnitude at r ~ 70AU, because of the reduced 
heating. For both the models FO.l and FC, EUV dominates for r < lOOAU, and beyond 
~ 100 AU FUV and dust collisions heat gas in the flow region. The increased heating causes 
a factor of ~ 20 change in Mpe at r ~ 100 AU between Models FC and Model O. At the 
outer edge of the disk (r ~ 200 AU), although FUV and dust collisions contribute almost 
equally to the heating for Model FC, there is not much change in Mpe when compared to the 
purely optical case (Model O) with only dust heating of gas. The gas temperature at the 
flow surface is marginally higher for Model FC, but this is compensated by the increased gas 
density at the flow surface in Model O, resulting in only a factor of ~ 3 change in Mpe at 200 
AU. We therefore find that FUV photoevaporation can cause significant mass loss at early 
stages of high accretion and associated high FUV fluxes, but at later stages when accretion 
has ceased, chromospherically generated FUV is not very significant in dispersing the disk. 

Our steady state analysis of the disk around a IM0 star with standard X-ray luminosity 
suggests that FUV/X-ray photoevaporation may not create a gap in the disk as in EUV 
photoevaporation. Although the mass loss rate peaks at ~ 3 — 40 AU in Model S, we 
find that locally M^g is always lower than the corresponding Mace-, indicating that viscous 
replenishment is always faster and therefore that a gap may never form. For the same disk, 
but at different epochs with higher and lower FUV fields, the model results show no evident 
peak in the FUV-generated mass loss rate. Thus it appears that in the standard case, a 
gap almost appears when Mace ~ 3 x 10~^ M0yr~^, but not quite. However, note that 
our models are calculated for fixed epochs and with the standard X-ray luminosity. Time- 
dependent models that include viscosity and futher explore the range of X-ray and FUV 
luminosities are needed to fully assess the manner in which the disk is dispersed by EUV, 
FUV and X-ray disk photoevaporation. It may be possible that such models will reveal gap 
formation at ~ 10 AU from the combined effects of FUV/X-rays. In fact, we find below in 
§3.3 that enhanced X-ray luminosities might lead to such a gap. 
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3.3. X-rays 

We have up to now treated the combined effects of FUV and X-ray photons in heating 
the disk gas and causing neutral photoevaporative flows. We next examine in some detail 
the effect of X-rays on disk heating and photoevaporation. We consider cases where we vary 
the stellar X-ray luminosity by factors of 10 from the fiducial model. In order to isolate the 
effect of X-rays on disk evolution, we also consider a hypothetical situation where the star 
produces no EUV or FUV photons. 

Figure [5] shows the results for our three models where M^e is shown as a function of 
radius. Note that all three runs (S, XIO, XO.l) also include EUV and FUV at our fiducial 
values (Tabled]). Increasing the X-ray luminosity increases the gas temperature which in 
turn increases Mpe, as expected. The effects of X-rays are most pronounced in the inner 
disk (r < 30AU) where they directly contribute to gas heating. X-rays in combination with 
FUV can heat the gas to ~ 1000 — 5000K in the inner disk (r < 1 — 10 AU), and they can 
significantly increase the mass loss rates here, by a factor of ~ 3 for a factor of 10 change 
in X-ray luminosity. Note that FUV still dominates the heating here (~ 75%). However, 
at small disk radii (r < 3AU), EUV-initiated photoevaporation is higher than that due to 
FUV and X-rays. The radius at which X-ray heating becomes important decreases as the 
X-ray luminosity increases, and varies from 1 — 5 AU for the range of X-ray luminosities we 
consider here. X-rays also affect the photoevaporation rates in the outer disk, by a factor 
of ~ 3 when we increase the X-ray luminosity by a factor of 10 (Model XIO) from the 
standard case (Model S). Decreasing the X-ray luminosity by a factor of 10 (Model XO.l) 
does not affect Mpe in the outer disk significantly, as gas heating here is mainly due to FUV 
photons. In the inner r < 10 AU disk, however, decreasing the X-ray luminosity has a 
more pronounced effect on lowering Mpg, which will lower the possibility of gap creation by 
FUV/X-ray photoevaporation. 

We find that stellar X-rays are largely ineffective in directly dispersing the disk (also see 
Alexander et al. 2004), although they can significantly affect the heating and chemistry in 
the gas. In our test run with only X-rays but no FUV or EUV (Model X) , X-rays do heat the 
surface of the gas to ~ few lOOOK, but the low densities in these surface regions of the disk 
result in low mass loss rates from the disk (dashed line in Fig. [5]). At r ~ 1 — lOAU, where X- 
rays were earlier found to significantly affect Mpe in our disk models XIO, S and XO.l, we find 
that removing EUV and FUV irradiation has a very dramatic effect. Our model disks S and 
X have the same X-ray luminosity, but zeroing the UV flux reduces Mpe by a factor of ~ 100. 
We therefore conclude that X-rays by themselves do not cause significant photoevaporation. 
We find that X-rays contribute to the photoevaporation in an indirect manner, due to the 
fact that X-ray ionization contributes electrons which thereby increase FUV heating. FUV 
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heating is mainly by the grain/PAH photoelectric heating which increases as the gas electron 
abundance increases, since the electrons reduce the positive charge of the grains/PAHs (e.g., 
Bakes & Tielens 1994, Weingartner & Draine 2001). Therefore X-rays act to amplify FUV 
heating by increasing the electron abundance in the gas. Even without FUV, Fig. O in case 
X shows a sudden rapid rise in Mpe at r ~ 75 AU, where Zfiow shifts from low density (n ~ 
lO^cm"^) X-ray heated gas to higher density {n ~ lO^cm"^) gas heated primarily by collisions 
with dust (~ 70% of total heating, 30% by X-rays). This is clearly seen by comparison with 
the model disk O (with only optical photon heating) and the addition of X-rays (Model X) 
increases mass loss only by a factor of ~ 2 due to the small increase in temperature. Stellar 
optical photons, via dust absorption and collisional heating of gas, are primarily responsible 
for the rise in Mpe at r > 75 AU in Model disk X, and X-rays make only a small contribution 
to the mass loss here. The main effect of X-rays on disk photoevaporation is thus indirect, 
with little direct photoevaporation caused by X-ray heating of gas. 

Although ineffective in removing the bulk of the mass in disks. X-rays may have a 
significant effect on its evolution. In our disk model with increased X-rays (Model XIO) 
Mpe ~ Mace in the inner r ~ 10 AU region of the disk, with Mpe ~ 2 x lO'^Mgyr-^ This 
suggests that a combination of FUV with high X-rays may be successful in driving a gap 
in the inner region. However, as discussed earlier, time-dependent models are needed to 
validate this tentative result. 

An earlier study by Alexander et al. (2004) found X-ray photoevaporation to be unim- 
portant. However, they did not include the effects of FUV heating in their models and only 
consider direct heating by X-rays and the resulting photoevaporation. They find that the 
upper limit to the rate of change of surface density by X-rays is ~ 10~^^ g cm~^ s~^ at 10 
AU, comparable to the EUV rate at that radius, and that most of the X-ray flux is absorbed 
close to the inner edge of the disk. A more recent study by Ercolano et al. (2008) suggests 
that X-ray photoevaporation may in fact be quite significant, with mass loss rates ~ 10~^^ g 
cm~^ s~^ at 10 AU. However, they ignore the effects of FUV and EUV, and do not accurately 
calculate the disk vertical structure. Note that the mass loss rates are very sensitive to the 
density at the base of the flow (which falls rapidly with z as e~^^^'^^^ , H being the disk scale 
height) and to the gas temperature (exponential dependence of S, Eqs. ISM]), necessitating 
detailed and self-consistent disk structure models. Our results are more in agreement with 
the results of Alexander et al. (2004), and we find that X-ray photoevaporation in itself is 
not very effective in dispersing the disk. We emphasize that X-rays can nevertheless affect 
the mass loss rates from the disk via the above indirect effect of raising the degree of ion- 
ization in the disk and increasing the efficiency of FUV-induced grain photoelectric heating, 
and in cases with high X-ray luminosities, promote the formation of gaps in the inner r ~ 10 
AU disk. 
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3.4. Dust Opacity and Surface Density Distribution 

Dust opacity is an important parameter influencing disk structure, and we investigate 
how tliis affects disk photoevaporation rates. As disks evolve, dust grains are believed to 
coagulate and grow and perhaps settle (e.g., van Boekel et al. 2005, Natta et al. 2007) 
reducing the opacity per H nucleus in the disk. FUV photoevaporation originates close 
to the Ay = 0.3 layer, and reducing the dust opacity may be expected to increase the 
photoevaporation rate due to the associated increase in gas density at the Ay = 0.3 layer. 
On the other hand, lower dust opacity implies lower FUV heating by the grain photoelectric 
effect, which would decrease the photoevaporation rate. Dust settling also causes a reduced 
flaring of the disk (Dullemond & Dominik 2005), and hence decreases the amount of high 
energy radiation intercepted by the disk, perhaps decreasing the photoevaporation rate. We 
study these effects by varying the dust cross section per H atom. We keep our dust grain size 
distribution flxed (from 50 A-20/im), and model a reduction in dust opacity by increasing the 
gas/dust mass ratio by a factor of 100. We also decrease the PAH abundance in proportion 
to the reduction in dust opacity. We flnd that r^isk decreases with a reduction in dust opacity 
(Fig. [6]), suggesting that as the disk evolves and grains settle and grow, the disk might lose 
its gas even more rapidly. However, the effect is not large, as seen in Fig. [6], with r^isk only 
increasing by a factor of 2, for a factor of 100 increase in dust opacity. Note that settling 
timescales depend on disk radius, and are typically longer in the outer disk where FUV 
photoevaporation dominates. Our derived disk lifetimes are typically short, ~ 1 Myrs, so 
that dust in the outer disk may not have adequately settled or undergone grain growth to 
signiflcantly affect disk survival times. 

We also consider a model with a steeper surface density distribution proflle with radius 
(S(r) oc r~^) compared with the standard disk model (S(r) oc r~^). We flnd that a steeper 
density distribution increases the disk lifetime, and that this effect is due to the resulting 
change in the disk structure. Steeper surface density distributions cause less flaring of the 
optically heated dust disk (e.g., Dullemond 2002), which results in a smaller angle being 
subtended by the disk and thereby reduced interception of the FUV and X-ray flux from 
the star. We note that even steeper proflles may result in the outer disk being partially or 
completely shadowed by the inner disk, perhaps limiting the rate at which mass loss can 
occur. Photoevaporation of the outer disk may occur until the surface density distribution 
is steeper than at which self-limiting point the outer disk is shielded by the inner disk, 
and the photoevaporation rate signiflcantly declines. However, in our present 1-l-lD dust 
radiative transfer model, we cannot investigate this effect and future 2D models are needed. 
Disks may also be shielded from incident FUV radiation if there are corrugations or ripples 
in the density distribution caused, for example, by spiral density waves. We do not consider 
such effects at present. We also note that future time-dependent models are needed to 
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investigate how the surface density distribution evolves as it loses mass due to FUV/X-ray 
photoevaporation, and how Mp^ and Tdisk are affected by the changing surface density profile 
in the disk. 



3.5. Disk lifetimes and stellar mass 

We next investigate the possibility that the disk lifetime is dependent on the mass of 
the central star. If disk mass loss rates are determined by FUV photoevaporation due to 
the central star, a peak in disk survival time with mass could possibly arise. Low mass stars 
may have lower initial disk masses than high mass stars, and their disks are more weakly 
bound by the stellar gravitational field. Disk lifetimes may therefore increase with stellar 
mass for low to intermediate mass (M* < 3Mq) stars. For massive stars (M,,, > TM©), most 
of the photospheric stellar emission is in the FUV ( 7 <M=i, < 2OM0) or in the EUV 
(M^, > 25M0) and their FUV or EUV luminosities are high. This may lead to very short disk 
lifetimes due to photoevaporation. It may be expected that disk lifetimes gradually increase 
with M^, peak at intermediate stellar masses (~ 3M0) and then rapidly drop at higher 
masses (M,, > TM©). Disk lifetimes may have a direct bearing on the likelihood of planet 
formation in the disk, which may then be a function of stellar mass. We test the hypothesis 
that disk lifetimes show a peak with M^, with models of disks around stars of different stellar 
masses, and show below that although disks around high mass stars are short-lived, those 
around stars from 0.3 — 3M0 have approximately equal lifetimes. 

The FUV luminosity of stars is a highly variable quantity for stars of low mass where 
the FUV does not arise from the photosphere, and can range over orders of magnitude 
for a given stellar mass. In lieu of any definitive observational data characterizing FUV 
luminosity as a function of stellar spectral type or mass, we use a mass-accretion defined 
FUV flux as described in GH08, i.e. we use the empirically derived accretion rate vs. mass 
relation of MuzeroUe et al. (2003, 2005) and assume that Mace ^ M^. We normalize this 
to the "median" accretion rate for a 1 Myr old solar mass star. Mace = 3 x 10^® Mq/jt. 
We then convert Maec into an accretion luminosity following the procedures of Gullbring 
et al. (1998). For massive stars, a significant portion of the FUV flux arises from the 
photosphere. It should be noted that due to the uncertain nature of the accretion-generated 
FUV spectrum, our results are inherently of a qualitative nature. Table [2] summarizes all 
the parameters we adopt for the different stellar masses. We additionally include a 3OM0 
case where the disk photoevaporation rate is entirely driven by the EUV radiation from 
the star and use the analytical results of HJLS94 to estimate disk lifetimes. Recall that 
we assume the initial disk mass Mdisk is proportional to the stellar mass, Mdisk = 0.03M^,. 
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This assumption together with the empirical criterion Mace oc M^, imphes that the viscosity 
parameter a increases with stellar mass (see Table 2). 

Figure [7] shows the calculated mass loss rates Mp(.{r) for the range of stellar masses 
considered here. Evidently, Mpf,{r) does not monotonically increase with stellar mass, even 
though the FUV luminosity is an increasing function of stellar mass. There is a peak at the 
gravitational radius in each case, which shifts to larger radii for more massive stars. The 
gravitational radius Vg increases with the mass of the star. This implies that in order to 
escape the higher gravitational fields, gas temperature in the flow needs to be higher at a 
given radius for higher mass stars. Gas temperature, however, is only a weak function of 
the strength of the FUV field. Therefore, lower mass stars have lower gravitational fields 
that only loosely binds their lower disk mass, but even with their lower FUV luminosities 
can heat the gas sufficiently to cause significant photoevaporation. On the other hand, 
the luminosities of very massive stars are photospheric in origin and significantly higher 
than accretion-generated FUV luminosities of intermediate and low mass stars, which are a 
fraction of the bolometric luminosity. Strong EUV and FUV fiuxes from massive stars can 
therefore erode their disks rapidly. 

We find that disk lifetimes Tdisk do not depend on the mass of the central star for 
0.3 — 3M0 stellar masses. We calculate Tdisk as described in §3.1 for each disk model. 
Figure [7] also shows the disk lifetime Tdisk as a function of stellar mass for our small sample 
of disk models. However, Tdisk drops sharply for both massive stars we consider, the 7Mq 
star and the 3OM0 star with EUV photoevaporation. These disks are very short-lived with 
Tdisk ~ 10^ years. Our result of a short disk survival time (~ 10^ years) for the most massive 
stars is supported by the observational evidence of the rarity of disks around stars more 
massive than ~ 7Mq (e.g., Fuente et al. 2007, Manoj et al. 2007). Massive stars may indeed 
lose their disks very rapidly due to their strong FUV and EUV fields. We do not find Tdisk 
increasing with stellar mass from 0.3 — 3M0 mainly because the FUV photons that dominate 
photoevaporation are generated by accretion, characterised by the relation Mace oc M^. Note 
that the heating at the fiow surface in the O.3M0 case (which has very low FUV) is by optical 
heating of dust grains which then collide with and heat gas. Therefore, even though low 
mass (M^ < IM0) stars start with less massive disks, and these disks are weakly bound, the 
much lower FUV for the lower mass stars also decreases photoevaporation rates and leads to 
lifetimes comparable to those for 1 — 3M0 stars. We conclude with the caveat that this result 
is very dependent on the various assumptions made in the analysis. Although our choice of 
model input parameters was guided by observational constraints, in many cases the existing 
data is not definitive. A different choice of input parameters, for example, initial disk mass 
proportional to M^, or Mace oc M*, may result in a peak in disk lifetimes with stellar mass. 
Our analysis involves assumptions that affect the value of a in the disk, affecting tjy, and 
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hence the derived disk lifetimes. Furthermore, a may vary both as function of disk radius 
(r) as well as disk height (z) due to varying ionization levels caused by MRI activity, and 
we consider a constant a throughout the disk. Our results are hence qualitative in nature. 

Observational evidence for a peak in disk lifetimes as a function of mass is not very 
conclusive. There appears to be a peak in disk fraction as a function of stellar type (a 
surrogate for stellar mass) in some clusters (Lada et al. 2006; Luhman ct al. 2008; Strom, 
private communication), as well as an increased incidence of debris disks with increasing 
stellar mass (Trilling et al. 2007, Hillenbrand et al. 2008, Cieza et al. 2008). Carpenter 
et al. (2006) report disk fractions as a function of stellar mass for the 5 Myr Upper Sco 
association and they find that that at least the inner disks (as probed by mid-infrared dust 
emission) are longer lived for low mass stars (0.1 — 1.2 Mq), and for more massive stars 
(1.8 — 20 M0) and are shortest in the intermediate 1 — 2 Mq range. On the other hand, 
observations of disks in young clusters such as Tr 37 and rj Cha (Sicilia-Aguilar et al. 2006, 
Megeath et al. 2005) indicate that disk fractions increase monotonically with central star 
mass (K-M spectral type range), and are less likely to be found around lower mass stars. 
Recent studies of associations in Orion (Hernandez et al. 2007) indicate a peak at a spectral 
type of MO, which corresponds to a mass of O.SM©. Damjanov et al. (2007) conclude that 
the disk frequency in Chameleon I is not a function of stellar mass, for a spectral type range 
from K3-M8. An important caveat is that the observational data report dust detections and 
do not directly trace the gas component or mass of the disk. Conversion of dust observations 
to disk masses make assumptions about the dust opacity law which is rather poorly known 
(e.g., Andrews & Williams 2007) and also assume a gas/dust mass ratio. 

Several studies may help validate or rule out photoevaporation as the principal agent 
in disk dispersal. Further observations of dust disk fractions in clusters may help resolve 
the seemingly contradictory results obtained so far. New studies that probe gas directly 
will be helpful. Better theoretical models are also warranted. The inclusion of viscosity in 
time-dependent disk models (Gorti, Dullemond & HoUenbach, in preparation) and perhaps 
improved determinations of UV fluxes to input into these models may also help reconcile the 
present differences between observational data and model predictions. 

4. Summary and Conclusions 

We present results of a study on the effects of optical, FUV, EUV and X-ray photons 
from a star illuminating its own circumstellar disk and their role in dispersing the disk gas. 
Earlier work on photoevaporation focused on the ionizing EUV photons either from a nearby 
massive star or from the central star in an early energetic phase in stellar evolution. Here we 
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mainly examine the role of FUV and X-ray photons which are produced even by low-mass 
stars during and after their accretion phase and their effectiveness in driving mass loss from 
the disk. We include optical, EUV, FUV and X-ray heating from the central star and find 
that FUV photoevaporation is most efficient in getting rid of bulk of the disk gas, located 
in the outer regions of the disk. 

We find that FUV photoevaporation induces neutral flows in the disk that are most 
significant at the outer disk regions (r > lOOAU), where most of the disk mass resides. 
Typical mass loss rates for a IMq star are ~ 3 x 10^^ M0yr^^ at 100 — 200 AU, and FUV 
photons incident on the disk (of initial mass 0.03M*) can therefore deplete most of its mass 
rapidly, on timescales of the order ~ 10^ yrs. X-rays do not cause significant flows, but can 
enhance FUV photoevaporation rates by ionizing the gas and increasing the efficacy of FUV 
grain photoelectric heating of gas. EUV photoevaporation is effective in the inner (r < 3AU 
for a IM0 star) regions of the disk (Mpe ~ 10"'^'^ Moyr"^) where the stellar gravitational 
potential is high, and higher temperatures (10'* K) are needed for escaping flows. EUV 
photoevaporation may also affect disk evolution at later stages by creating gaps in the disk 
(at r ~ 7 AU, IM0 case). FUV radiation acting with relatively high X-rays could potentially 
create gaps at 3 — 30 AU at some stages of evolution, but this is not very conclusive from 
the static model results of this paper. 

FUV photoevaporation depletes the disk mass early in its evolution, before EUV photons 
arc able to penetrate the protostellar winds that accompany accretion. Wc suggest that disk 
dispersal is initiated by FUV photons (aided by X-rays) which truncate the disk at the radius 
where viscous timescales equal photoevaporation timescales (typically at r ~ 100 AU) and 
remove much of the disk mass. However, time dependent models that follow the decrease 
in FUV luminosity as the accretion rate drops in time are needed. As the disk mass drops, 
viscous accretion rates decrease to values where EUV photoevaporation can create a gap. 
This occurs when the EUV-induced mass loss rate at Vgji ~ 7(M^/(1Mq) is higher than 
the accretion rate here. The inner disk then drains viscously and the direct illumination 
of the inner rim (at ~ 7AU) by EUV, FUV and X-rays then results in photoevaporation 
of the remaining outer disk gas. However, it is the FUV photoevaporation of the outer 
disk mass reservoir that sets the disk lifetime. We flnd that FUV photoevaporation rates 
(and hence disk hfetimes) are sensitive to the FUV luminosity of the star. If the FUV 
luminosity is mainly generated by accretion activity, the mass loss rates decline with time 
as accretion diminishes in the disk. FUV generated by stellar chromospheric activity is not 
very significant in dispersing the disk, and produces mass loss rates similar to that by stellar 
optical photons and dust heating of gas at ~ 200 AU. Disk photoevaporation timescales 
{Tdisk) are calculated at ~ 1 Myrs after the penetration of the disk wind by FUV and X-rays 
(at an age ~ IMyr) and we estimate that disks may last ~ 2 Myrs. 
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We also find that the disk dust opacity may influence the photoevaporation rates by 
affecting gas density at the flow origin and the degree of flaring. For a factor of ~ 100 
reduction in dust opacity, we flnd that Tdisk decreases by a factor of ~ 2. The surface density 
distribution of gas in the disk may also affect the photoevaporation rate, as steeper radial 
distributions decrease the radiation intercepted by the disk and reduce the mass loss rates; 
changing our surface density from S oc r^^ to r^^ increases the disk lifetime by ~ 50%. 

Disk lifetimes do not appear to depend on the mass of the central star for stellar masses 
between 0.3 and 3Mq stars. Here, we derive lifetimes of ~ 10^ years. High mass stars 
(M* > 7M0) lose their disks rapidly in ~ 10^ years because of their very high photospheric 
FUV and EUV fields. 

Our derived disk lifetimes of ~ 10^ years after the FUV penetration of the protostellar 
wind for typical solar-mass stars are consistent with observations, although this study does 
neglect dynamics and only crudely considers viscosity. In addition, derived disk lifetimes 
may be longer if the initial disk mass is greater than 0.03M*. However, our qualitative 
estimates of the "mean" disk lifetime by assuming "mean" or average stellar properties at a 
typical age of 1 Myr, compare very well with observationally inferred disk lifetimes, without 
requiring very high EUV photon fluxes (Alexander et al. 2006) or special birthplaces near 
massive OB stars in a cluster environment (e.g., HoUenbach et al. 2000, Adams et al. 2004, 
Johnstone et al. 2004). 

We offer a concluding important caveat. The present calculations are simplistic and 
only quahtatively include the effects of viscosity. As the FUV fleld, accretion rates, and disk 
properties evolve with time, a proper treatment needs to consider all these effects in a time- 
dependent model. A future goal will be to incorporate the present FUV, EUV and X-ray 
photoevaporation models into a self-consistent and time-dependent solution that includes 
viscosity (Gorti, DuUemond & HoUenbach, in preparation) 
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Table 1: Fiducial Disk Model - Input Parameters 



Disk mass 0.03 M* 

Surface density E(r) oc 

Inner disk radius 0.5 AU 

Outer disk radius 200 AU 

Gas/Dust Mass Ratio 100 

Dust size distribution n{a) oc a~^'^ 

Min. grain size Omm 50A 

Max. grain size Umax 20/xm 

(Th 2 X lO-^^cmVH 

PAH al)uii(laiice/H 8. 1 x IQ-'" 



Table 2: Standard Stellar Input Parameters as a Function of Mass 
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Log 0EUV 
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(erg s~"^) 


(erg -^) 




(Mo) 




0.3 


2.3 


3360 


0.55 


2.7 (-9) 


30.3 


29.6 


39.9 


0.009 


0.005 


0.5 


2.12 


3771 


0.93 


7.5 (-9) 


30.9 


29.8 


40.1 


0.015 


0.009 


0.7 


2.54 


4024 


1.72 


1.5 (-8) 


31.3 


30.2 


40.5 


0.021 


0.012 


1.0 


2.61 


4278 


2.34 


3.0 (-8) 


31.7 


30.4 


40.7 


0.03 


0.02 


1.7 


3.30 


4615 


5.00 


8.7 (-8) 


32.3 


30.7 


41.0 


0.051 


0.036 


3.0 


4.83 


5004 


14.85 


2.7 (-7) 


32.9 


28.7 


39.0 


0.09 


0.066 


7.0 


3.22 


20527 


1687 


1.5 (-6) 


36.5 


30.8 


44.1 


0.21 


0.07 



-31 - 




r (AU) r (AU) 



Fig. 1. — Parameters at the flow surface of a photoevaporating disk. In the left panel, the 
location of the ionization front and the neutral photoevaporative flow surface are marked. 
The base of the neutral FUV-flow lies slightly above the Ay = 1 layer throughout the disk. 
The right panel shows the temperature and density at the base of the flow. For r < 50 
AU, the gas temperature falls steeply with increasing r, whereas it is a shallow function of r 
beyond 50 AU. The gas density profile shows a sharp kink at r ~ 3AU, where the fiow shifts 
from the EUV to the FUV heated region, but is otherwise fairly shallow through the entire 
disk. 



-32- 





— R 


_i 


1 1 
Model 

IVL vJ U 1 


1 1 1 1 1 
S: Log Lpuv=31.7; 

VC- I rta T — ?Q Q- 


1 1 1 1 1 1 1 l_ 

Macc=3xl0-° yr-' 

IVl — \J 

acc _ 




-7 








Standard _ 
FUV/X-rays 


1 






Mace 




/S 










/ 


© 


-8 








V / cr* 
\ / ■ ^ 














cv? 
CV} 


-9 








\ / .'Chromospheric 
\ / /FUV/X-rays ~ 














0) 

o. 








EUV 




Log 


-10 
-11 


1 


1 . 


1 1 1 1 1 


EUV ^' _ 

1 1 1 1 1 1 1 1 - 











1 


2 



Log r(AU) 



Fig. 2. — The photoevaporation rate expressed as a local mass loss rate (27rr^S) and shown 
as a function of radius for our fiducial disk around a IMq star (Model S). Inside of 3 
AU, EUV photoevaporation dominates, and at r > 3 AU, FUV/X-ray heating dominates 
photoevaporation. The mass loss rate peaks at the outer edge of the disk, and shows a dip at 
around 40 AU, where the disk is expected to last longest. There is also a peak at ~ 10 AU, 
where a gap may be created at favourable epochs in disk evolution, when Mpg > Mace, the 
viscous accretion rate (dot-dashed fine). The figure also shows Mp^ for another disk model, 
where accretion has ceased and the FUV is generated by the stellar chromosphere (Model 
FC). In this case, EUV photoevaporation rates are higher than FUV for r < 90AU, beyond 
which FUV photoevaporation dominates. 
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Fig. 3. — The timescale for mass loss due to photoevaporation {tpe ~ E(r)/il(r)) for our 
standard disk as a function of radius (solid line). Also shown is the viscous timescale in the 
disk (dashed line). The disk lifetime r^isfc is estimated as the intersection of tp^ and at 
some truncation radius r^. In a time Tdisk, the disk mass interior to is drained by viscosity 
whereas the disk outside of rj photoevaporates. 
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Fig. 4. — The photoevaporation rates for models with accretion-genererated FUV luminosi- 
ties of 4 X 103°(Model FO.l), 4 x lO^i (Model S), and 4 x lO^^ erg s'^ (Model FIO), and 
corresponding mass accretion rates of 3 x 10~^,3 x 10~^ and 3 x 10~^ Mq/jt respectively. 
Also shown are a case with no accretion, but only the chromospheric FUV component of 
2 X 10^° erg s~^ (Model FC), and a model with only optical photons and no UV or X-rays 
(Model O). For models FO.l and FC, EUV photoevaporation dominates for r > 100 AU. 
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Fig. 5. — The photoevaporation rate for three different X-ray himinosities, 2 x 10^^ erg s ^ 
(Model XIO), 2 x lO^o erg s^^ (fiducial case, Model S) and 2 x fO^^ erg s^^ (Model XO.l) as 
solid lines. All these runs also include FUV (2 x 10^^ erg s~^) and EUV (2 x 10^° erg s~^) 
illumination of the disk. The dotted line shows a hypothetical case (Model X) with no EUV 
or FUV, but only an X-ray luminosity of 2 x 10^° erg s~^. The rapid rise in Mpe for model X 
is not caused by X-rays, but is mainly due to optical photons incident on dust grains which 
collide with and heat the gas at the disk surface. Disk model O with no FUV or X-rays is 
also shown for comparison. 
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Fig. 6. — Disk lifetime, Tdisk for the IM© case, with different dust opacities in the disk, and 
the surface density distribution S(r) oc r~^. Also shown is a case with a dust cross section 
of 10~^^ cm~^ per H nucleus and with a steeper surface density distribution, E(r) oc r~^. 
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Fig. 7. — Disk mass loss rates Mpg and derived disk lifetimes are shown for a range of stellar 
masses. Disk lifetimes are determined by finding the radius where the photoevaporation 
timescale is equal to the viscous timescale t,^, and adopting this timescale (see text). 
Disk lifetimes are nearly independent of stellar mass for < 3Mq and sharply decrease for 
> 3M0. 



